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Abstract 



r^ , Consider a random walk whose (light-tailed) increments have positive mean. Lower 

and upper bounds are provided for the expected maximal value of the random walk 
until it experiences a given drawdown d. These bounds, related to the Calmar ratio 
in Finance, are of the form {exp{ad} — l)/a and (if exp{ad} — l)/a for some K > 1, 
in terms of the adjustment coefficient a {E[exp{—aX}] = 1) of the insurance risk 
literature. Its inverse — has been recently derived by Aumann & Serrano as an index 
of riskiness of the random variable X. 

This article also complements the Lundberg exponential stochastic upper bound and 
\^ , the Cramer-Lundberg approximation for the expected minimum of the random walk, 

CO ' with an exponential stochastic lower bound. The tail probability bounds are of the form 

r — I . C exp{—ax} and exp{—ax} respectively, for some -^ < C < 1. 

f^ ' Our treatment of the problem involves Skorokhod embeddings of random walks in 

OO , Martingales, especially via the Azema-Yor and Dubins stopping times, adapted from 

^^ ' standard Brownian Motion to exponential Martingales. 

AMS classification: Primary 60G50 , 60G44 ; secondary 91B30 
S^ . Keywords and phrases: Calmar ratio, Cramer Lundberg, Drawdown, Random Walk, 

H ' Skorokhod embeddings 

1 Introduction 

Drawdowns of Brownian Motion with positive drift. Let {W{t) | t > 0, W{0) = 0} 
be Standard Brownian Motion (SBM) and let {B{t) \ B{t) = fit + aW{t), t > 0} be 
Brownian Motion (BM) with drift fj, > and diffusion parameter a £ (0,oo). For d > 0, 
define the stopping time 

rf ^ = min{t| max B{s) > B{t) + d} (1) 

0<s<t 

to be the first time to achieve a drawdown of size d. That is, t^ is the first time that 
BM has gone down by d from its record high value so far. As motivated by Taylor [20], 



an investor that owns a share whose value at time t is V^ = VQexp{B{t)), may consider 
selhng it at time t^^ (for some d > 0) because it has lost for the first time some fixed 
fraction 1 — exp(— d) of its previously held highest value Vq exp{M(i) (where M^ = Mf^ = 
maxQ^^^^BM B{s) = B{t^^^) + d), a possible indication of change of drift. 

As pointed out in Meilijson [17], drawdowns are gaps for Dubins &: Schwarz [12]), extents 
for Goldhirsch &: Noskovicz [13]) and downfalls for Douady, Shiryaev & Yor [9]). Taylor 
[20] (see also [17]) presents a closed form formula for the joint moment generating function 
of T^^ and B{t^^), from which it follows that M^^ is exponentially distributed, with 
expectation 

E[Mi^] = ^{eM^^d}-l). (2) 

Maximum of Brownian Motion with negative drift. The maximum 

max(i?M) = mit>o{B{t)} is well known to have the exponential distribution 

P(max(BM) > x) = 1 A exp{ ^x}. (3) 

This article contributes to the generalization of (2) and (3) from BM to random walks 
(RW). There is a rather vast literature on the maximum of RW with negative drift. Kingman 
[15] showed that P{uiayi{RW) > x) ~ 1 A exp{ — -^x} for small /i, Siegmund [18] studied 
first order corrections to this approximation via renewal-type overflow distributions and 
Chang Sz Peres [8] developed asymptotic expansions of P{inax{RW) > x) for the Gaussian 
case. Blanchet &: Glynn [6] improved on these approximations. In the insurance risk 
literature, exponential bounds and approximations of P{niax{RW) > x) are referred to as 
Lundberg's inequality or Cramer-Lundberg approximations (see Asmussen's comprehensive 
treatise [1]). 

This paper is methodologically different from the above; instead of relying on change 
of measure and renewal theory, our setup involves exponential martingales and Skorokhod 
embeddings, in a way reminiscent of Wald's [21] method for deriving the OC characteristic 
of the Sequential Probability Ratio Test. As part of the change, we will give up on trying to 
save the inaccurate role of -^ as the exponential rate in the questions under study, in favor 
of the so-called adjustment coefficient of the insurance risk literature, provided by the a 
solving i?[exp{— aX}] = 1. However, the rate -^ will stay around: the RW will be coupled 
with a BM for which -^ is a. 

More explicitely, using Skorokhod ([19] and also [4, 7, 10, 16]) embeddings, mean-zero 
RW can be viewed as optional sampling of SBM. This idea will be mimicked here to embed 



the exponential Martingale exp{— qS'„} into the Martingale exp{—aB{t)}. This method 
could be useful in obtaining other approximate extensions of pricing under log-normal 
models to more general distributions. 

Aumann &: Serrano [3] asked a scalar index of riskiness Q{X) of the random variable 
(r.v.) X to satisfy an homogeneity axiom Q{tX) = tQ{X) and a duality axiom that models 
the increased preference of a more risk averse individual for constant wealth w over random 
wealth w + X. The unique solution (up to a multiplicative constant) is the inverse — of the 
adjustment coefficient. The role played by a in our subject matter is clearly consistent with 
riskiness - a large a corresponds to low risk, as it (i) protects against heavy initial losses 
before eventual divergence of the RW to oo, and (ii) makes the RW reach high yield before 
experiencing sizable drawdowns. 

The Calniar ratio (see Atiya &: Magdon-Ismail [2] and the implementation of their work 
in the Matlab financial toolbox) of a financial asset with positive drift is a measure of the 
likely drawdown in the logarithm of its price in a given interval of time, such as a year. Since 
height (and time) are exponential in the drawdown, the Calmar ratio is heavily infiuenced 
by the length of this time interval. Besides, typical drawdown in a given time span is harder 
to analyze than our subject matter, typical height (or time) to achieve a given drawdown. 
We propose the use of the adjustment coefficient or its inverse as a Calmar-type measure of 
the risk of a financial asset, and provide simple approximate formulas to quantify its effects. 
The more commonly used Sharpe index, or ratio of net drift (drift minus market interest 
rate) to volatility (standard deviation), lets volatility penalize the asset even when it favors 
gains. In contrast, drawdown-based indices measure risk in a more reasonable asymmetric 
sense. 

2 Results 

From now on, we only consider BM and RW with positive drift and thus unify the presentation 
of the two problems, by switching from the commonly studied maximum of BM and RW 
with negative drift to the equivalent treatment of the minimum of BM and RW with positive 
drift. 

Consider a r. v. X r^ F with < E[X] < oo and P{X < 0) > 0. Assume further the 
existence of a > such that £^[exp{— qAT}] = 1. Since (if finite) the moment generating 
function *(t) = E[exp{tX}] is strictly convex with ^'(0) = E[X] > and ^(t) -^ oo 



as |t| -^ cxD, such a exists and is unique as long as the moment generating function is 
finite wherever relevant. This assumption is satisfied e.g. for Gaussian r.v.'s and for r.v.'s 
bounded from below. If X ~ N{fj,,a^), then a is indeed -^ (see (2)). 
Besides a, we need other characteristics of the distribution F. 



d+ = - sup -log(S[e-"(^-^)|X >x]) 

a 0<a:<eSF 



d- -- 


= - sup log(E[e"(^-^)|X<x]) 


do = 


= d-^ + d" 



(4) 



where esi;' = sup{y\F{y) < 1} and eip = mi{y\F{y) > 0} are the essential supremum 
and infimum of F. By Jensen's inequality, d'^ is bounded from above by the simpler and 
more natural sup^E[X — x\X > x] and d^ is accordingly bounded from below. These 
constants are defined in terms of excesses of the r. v. X itself, unlike the Siegmund or 
Cramer-Lundberg approximations, built in terms of the renewal overflow distribution of 
the random walk with X-increments. 

Let Xi ; i = 1,2, ... be i.i.d. F— distributed random variables and let 5*0 = ; Sn = 
^^=1 Xi be the corresponding random walk (RW). The definition of a makes exp{— a;S'„} 
a mean-1 Martingale. 

Drawdowns, maximal heights M^ achieved prior to drawdowns and the corresponding 
stopping times t^ can be defined for random walk in much the same way they are defined 
for Brownian Motion. Re-stating Lundberg's inequality as the RHS of (6), the purpose of 
this paper is to prove the other three inequalities in (5) and (6) 
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<E[Mr]< (5) 



a a 



For X > 0, 



^-aix+d ) < p^_ jnin(i?VF) >x)< e""^ (6) 

We thus have lower and upper bounds for £'[MJ^ ] whose ratio stays bounded as d 
increases, provided do is finite. These bounds clearly show that drawdowns are logarithmic 
in the highest value achieved so far, and precisely identify the exponential rate a at which 
the latter grows as a function of the former. Besides claiming bounds on the mean, the 



upper bound in (5) is in fact a stochastic inequality: M^^ is stochastically smaller than 
the exponentially distributed random variable M^^ ■ We do not provide a stochastic lower 
bound. In contrast, (6) provides stochastic upper and lower bounds on the minimum of 
RW. 

The upper bound in (5) can be improved by letting do depend on d and be defined as do in 
(4) but restricting the maximization to 

X £ (0, minxes F,d)) and x G {m.ax{eip,—d),0). This improved upper bound is finite for 
every d. 

3 A few examples 

Example 1: The Gaussian case. Let X have a normal distribution with positive 
mean /j, and standard deviation a. Let and <& stand respectively for the standard normal 
density and cumulative distribution function. Then 

and this proves the two following rather elegant formulas. 

%^ < E[M,] < ^--^^ (8) 

^(t^{^ ~^ - ^(-min(i?VF) > a;) < e"^"" (9) 

If we view the normal random walk as sampling Brownian Motion with drift fi and 
diffusion coefficient a at regular intervals, a is independent of the grid length S but the 
three d's are not, predictably vanishing with the grid length. 

The LHS of (7) is well known and easy to obtain from the formula exp{/i + cj^t^/2} of 
the moment generating function of the normal distribution. As for the RHS, it requires 
evaluating via 

^ ' ^ V2^ l-$(z) l-$(z) ^ ^ 

the expressions 

2 

E[e-''^\X>x] = E[e-'''■^'+''^^\n + aZ>x]=e~^^E[e^^^^\Z>'^^^] 
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(11) 



S[e-"-|X<x] = ^^ (12) 



from which the RHS of (7) follows, at least in the sense of plugging x = 0. To see that x = 
is indeed the correct choice for each side, observe that the normal distribution is IFR - has 
increasing failure rate (Mills' ratio i_^/ ) )■ But IFR implies that the residual distributions 
C{X — x\X > x) are ordered by stochastic inequality. Hence, so are the expectations of 
monotone functions, such as the exponential function. This argument applies equally to 
the two tails. 

Example 2: The double exponential case. Let X have density p6 exp{— Ox} for 
X > and (1 — p) fi ex.p{fix} for x < 0, with — ^ < p < 1. Then 

E[X] = ^-^^^; a=pf^-{l-p)e (13) 

with the corresponding bound ingredients 

Q ' (Ai + e)(i-p)' {\-p)e ^ ' 

The rate a exceeds y ^ J for all p ii fi > 0, but li n < 9 the opposite inequality holds 
for all p close enough to 1. 

Example 3: The shifted exponential case. Let X have exponential distribution 
with mean ^ shifted down by A < ^ so as to allow negative values and still preserve positive 
mean. It is easier to express the inverse function to a: 

-log(l + ^) = A (15) 

from which 

, n- - 1 - p-(^+a)A ^ 
(f+ = A • e°'^ = 1 + — ■ e°'^ = - • e""^" = H - flG) 



The Gaussian-motivated rate -,. ^ ^ is 29{1 — 9A), always smaller than a. That is, 
a random walk with shifted exponential increments gets to higher heights before a given 
drawdown than a normal one with the same mean and variance. 

Example 4: A dichotomous case. Let P{X = -1) = 1 — p and 

P{X = 1) = p > 2- Then a = log j^ and, obviously, d+ = d^ = 1. As is well known 
from the Gambler's ruin problem, the probability of reaching +1 before (integer) —d is 
^_^^^_~/'^,^\\ . M^ is nothing but the number of independent such attempts until a first 
"failure". Hence, it is (—1 plus) a geometric r. v., and its mean is 

E(M?^) = -1 + i = ^^(e"'^-l). (17) 

"- d J -r pj-obability 2p-r ' ^ ' 



For non-integer d, the ceiling of d should be substituted in (17). Even without doing so, the 
LHS of (5) is verified, because a = log j^ > ^^. To ascertain the RHS, take E{M^^) as 
worst-case ceiling and check that 2^^{e'^ — 1) is below the bound {exp{a{d+2)} — l)/a. 
Just as in the shifted exponential case, the rate of growth a exceeds the rate y ^ ' = 
2(Y:r^ — -) that would have obtained in the Gaussian case. However, 

Example 5: A skew dichotomous case. Let P{X = —1) = ^^ and 

P{X = 6(1 + e)) = j^. The mean is E[X] = ^ so let us take 
e = 0.2 to achieve positive mean and 6 = 0.1 to tilt the distribution towards bigger losses. 

Then a = 0.318 but -,. r , = 0.351. This shows that even for dichotomous variables the 

Var[X] 

inequality between the two can go both ways. Examples 2 and 5 suggest that yield-to- 
drawdown performance worse than Gaussian is obtained when the left tail is heavier than 
the right tail. 

In all the previous examples, the distribution F has non-decreasing failure rate and the 
"excess lifetime" over x looks shorter as x increases. That's why the d's are attained at 
X = (see (4)). It is easy to produce a four-point distribution with one negative atom in 
which d^ will be the distance between the two rightmost atoms. 

Example 6: A power-law right tail. If F is light left tailed but behaves hke power 

law at the right tail, then a is finite but d'^ is infinite because its maximand behaves like 

logx. To wit, 

E[e-»(^-^)\X >x] = ^ r -^^-—dt « ^ (18) 

^ ' ^ xJo (l + |)T+i ax ^ ' 

so — log(£'[exp{— a(X — 2;)}|X > x]) = log(x) -|- o(l). Although much smaller than 

E[X — x\X > x] = 0{x) (see the sentence following (4) and the remark in the next section), 

it still goes to cx). However, the improved definition of d^ sets it as log(d) up to a vanishing 

term. 

This example illustrates that yield-to-drawdown, while at least as high as the Brownian 

lower bound, may in principle be super exponential. 

4 Miscellaneous 

The record high value M^^ is exponentially distributed. This is so because as long 
as first hitting times of positive heights occur before achieving a drawdown of d, these times 
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are renewal times: knowing that M^^^ > x is the same as knowing that B has not achieved 
a drawdown of d by the time it first reached height x. But then it starts anew the quest 
for a drawdown. 

A direct argument for (2). Since the mean-1 Martingale exp{— aS} stopped at t^^ is 
uniformly bounded, it is also uniformly integrable. Hence, 

1 = i5;[e-^(-i'"')] = E[e-^^^'']e^^ (19) 

Since M^^ is exponentially distributed, i?[exp{— aM^^}] = 

5 Skorokhod embeddings in Martingales 

The problem as posed and solved by Skorokhod in [19] is the following: given a distribution 
F of a r. V. Y with mean zero and finite variance, find a stopping time r in SBM W, with 
finite mean, for which W{t) is distributed F. The Chacon-Walsh [7] family of solutions is 
easiest to describe: Express Y as the limit of a Martingale Yn = E\Y\J^n] with dichotomous 
transitions (that is, the conditional distribution of l^+i given J^n is a.s. two- valued), and 
then 

progressively embed this Martingale in M/^ by a sequence of first exit times from open 
intervals. 

Dubins [10] was the first to build such a scheme, letting J^i decide whether Y > E[Y] 
or y < E\Y] by a first exit time of W starting at E[Y] from the open interval (£^[y|y < 
-©[y]], £'[y|y > £'[y]]). it then proceeds recursively. E.g., if the first step ended at £'[y|y > 
£'[y]] then the second step ends when W, re-starting at £'[y|y > £'[y]], first exits the open 
interval {E[Y\E[Y] <Y< E[Y\Y > E[Y]]],E[Y\Y > E[Y\Y > E[Y]]]). 

One of the analytically most elegant solutions to Skorokhod's problem is the Azema- 
Yor stopping time Tj^y (see Azema & Yor [4] and Meilijson [16]), defined in terms of 
Hp{x) = E[Y\Y > x] = J^ ydF{y)/{l — F{x—)), the upper barycenter function of F, as 

Tay = minjt max W{s) > HF{W{t))}. (20) 

Among all uniformly integrable cadlag Martingales with a given final or limiting distribution, 
SBM stopped at the Azema- Yor stopping time to embed this distribution is extremal, in 
the sense that it stochastically maximizes the maximum of the Martingale (see Dubins &: 



Gilat [11] and Azema & Yor [4]). That is, if T^y embeds F then M^^^^ is stochastically 
bigger than the maximum of any such Martingale. 

The connection of the Azema- Yor stopping time to the Chacon- Walsh family becomes 
apparent (see Meilijson [16]) if the r. v. Y has finite support {xi < ■ ■ ■ < x^}. In this case, 
let J-'n be the a-field generated by min(y, Xn+i), that is, let the atoms of Y be incorporated 
one at a time, in their natural order: the first stage decides whether Y = xi (by stopping 
there) or otherwise (by temporarily stopping at £^[y|y > xi]), etc. This is precisely the 
Azema- Yor stopping rule: stop as soon as a value of Y is reached after having visited the 
conditional expectation of Y from this value and up. 

Clearly, there is a mirror-image notion Tay- to Azema & Yor's stopping time that 
stochastically minimizes the minimum of the Martingale. Simply put, apply Tj^y to embed 
the distribution of —X in —W. 

The stopping time Tjjay to be applied here is a hybrid of the Dubins and Azema &: 
Yor stopping times. It starts as the Dubins stopping time by a first-exit time of SBM 
W from the interval {E[Y\Y < E[Y]],E[Y\Y > E[Y]]). If exit occurred at the top it 
proceeds by embedding the law £(y|y > £'[!"]) by Tay in the remainder SBM starting at 
£^[y|y > -E[y]]. if, on the other hand, exit occurred at the bottom it proceeds by embedding 
the law £(y|y < E[Y]) by Tay- in the remainder SBM starting at E[Y\Y < E[Y]]. 

Once a distribution F is embeddable in SBM W, so is the random walk with increments 
distributed F. Plainly, embed Xi at time ri, then use the same rule to embed X2 at time 
T2 in the SBM W'{t) = W{ti + t) — W{ti), etc. Skorokhod's original idea was to infer the 
Central Limit Theorem for ^ from the Law of Large Numbers for ^^=^ . This idea was 
extended by Holewijn & Meilijson [14] from random walks to Martingales with stationary 
ergodic increments, to obtain a simple proof of the Billingsley & Ibragimov [5] CLT. 

Just as a random walk Sn can be embedded in SBM, the exponential Martingale 
exp{—aSn} can be embedded in the continuous-time continuous Martingale exp{—aB(t)}, 
where the BM B has drift /i and diffusion coefficient a such that -^ = a. At the time the 
RW reaches drawdown at least d, BM has also gone down by at least d, but may have gone 
higher in the meantime. Thus, rf is a.s. smaller than t^ . Now we may compute, under 
the obvious definition of A > d a.s., 

f,E[r,^''] = E[Bir,^^)]=E[Mi^]-d 



< ^xE[Tf'^]=E[B{Tf'^)]=E[Mf'^\-E[^] (21) 

so E[M^^] > ^[Mf^]. We have proved the LHS inequahty in (5) by the method of 
Couphng. 

The RHS inequahty in (5) is proved by using the Dubins - Azema & Yor stopping time 
'Tday for the above embedding. If this embedding in exp{— aVF} ends up with W below 
(resp. above) some x > (non-positive), the underlying process could not have reached the 
exponential barycenter height above (below) the closest support point to the right (left) of 
X, because then the stopped value would have been from this support rightwards (leftwards). 
For the first increment of RW following the maximal (minimal) value the relevant x is 0, 
but for values embedded later the starting x is lower (higher). It should now be clear that 
the BM path can't reach as far up as MJ^ + d^ nor as far down as B{tIP^) — d~ before 
RW achieves drawdown d. Hence, BM can't reach drawdown d + do before RW reaches 
drawdown d, or tJ^^ > t^ a.s. The RHS inequality in (5) follows. 

This RHS inequality holds stochastically, since the cumulative maximum of BM exceeds 
the cumulative maximum of the embedded RW timewise. Thus, it holds a fortiori if the 
former is measured later than the latter. This argument fails for the LHS because then the 
latter is measured before the former. 

Finally, a non-trivial stochastic upper bound on M^ can be defined even if do = cxd, 
as announced in the Introduction. In Example 5, this corresponds to adding a little over 
logd to the exponent, or, the upper bound is roughly d times the lower bound. 
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